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Abstract 

We discuss the impact of recent experimental results on the determination 
of atmospheric neutrino oscillation parameters. We use all published results 
on atmospheric neutrinos, including the preliminary large statistics data of 
Super-Kamiokande. We re~analyze the data in terms of both 
and Vf. channels using new improved calculations of the atmospheric 

neutrino flux. We compare the sensitivity attained in atmospheric neu- 
trino experiments with those of accelerator and reactor neutrino oscillation 
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searches, including the recent Chooz experiment. We briefly comment on 
the imphcations of atmospheric neutrino data in relation to future searches 
for neutrino oscillations with long baselines, such as the K2K, MINOS, 
ICARUS and NOE experiments. 

14.60.Pq, 13.15.+g, 95.85.Ry 
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I. INTRODUCTION 



Atmospheric neutrinos are produced in cascades initiated by collisions of cosmic 
rays with the Earth's atmosphere. Some of the mesons produced in these cascades, 
mostly pions and kaons, decay into electron and muon neutrinos and anti-neutrinos. The 
predicted absolute fluxes of neutrinos produced by cosmic-ray interactions in the atmo- 
sphere are uncertain at the 20% level. The ratios of neutrinos of different flavour are 
however expected to be accurate to better than 5%. Since z/g is produced mainly from 
the decay chain it —>■ followed by yU — > ev^j_Vei one naively expects a 2 : 1 ratio of z/^ 
to z/g. In practice, however, the expected ratio of muon- like interactions to electron-like 
interactions in each experiment is more uncertain [0,0. 

Several experiments have observed atmospheric neutrino interactions. Two under- 
ground experiments, Kamiokande and 1MB [0, use water-Cerenkov detectors. These 
experiments have detected a ratio of z/^-induced events to z/g-induced events smaller than 
the expected one In particular Kamiokande has performed separate analyses for 
both sub-GeV neutrinos and multi-GeV neutrinos [^, which show the same deficit. 
Although some of the experiments, such as Frejus |^ and NUSEX 0, have not found 
evidence for this anomaly, and others, e.g. Soudan2, are not yet conclusive, the recent 
Super-Kamiokande data provides a strong support for an atmospheric muon neutrino 
deficit. This encourages us to reconsider the analysis of atmospheric neutrino data from 
the point of view of neutrino oscillation interpretation. The recent improved data sample 
of Super-Kamiokande has now better statistics than achieved in the whole Kamiokande 
phase. 

We also include the new data of Soudan2 in our analysis. These new data as well 
as the previous experimental results are summarized in table |. Here i?^/e/-R^/^ denotes 
the double ratio of experimental-to-expected ratio of muon-like to electron-like events. 
The expected ratio -R*/^ is obtained by folding a prediction for the atmospheric neutrino 
flux with the properties of every individual detector through a Monte Carlo procedure. 

Apart from studying the impact of the new data, our motivation for the present 
reanalysis of atmospheric neutrino data is theoretical. In this regard, we first of all 
include the results of a recent calculation of the atmospheric neutrino fluxes as a function 
of zenith angle |jTO[, including the muon polarization effect Moreover, we develop 

an independent procedure for the comparison of results from different experiments. We 
demonstrate that our theoretical calculation of the energy distribution of the event rates 
is in good agreement with the MC expectations. The comparison of the experimental 
results presented below thus reflects the significance of the atmospheric neutrino anomaly 
and provides evidence for neutrino oscillations. 
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In this paper we analyze the impact of recent experimental results on atmospheric 
neutrinos from Super-Kamiokande and Soudan2 on the determinations of atmospheric 
neutrino oscillation parameters, both for the z/^ Vr and z/^ — > Vf. channels. In so-doing 
we take into account recent theoretical improvements in flux calculations and neutrino- 
nucleon cross sections. The new Super-Kamiokande data produce a downwards shift in 
the allowed (sin^ 26, Am^) region, when compared with pre-Super-Kamiokande results. 
Nevertheless we show that the t'^j oscillation hypothesis is barely consistent with the 
recent negative result of the Chooz reactor WM. The sensitivity attained in atmospheric 
neutrino observations in the z/^ — > channel is also compared with those of accelera- 
tor neutrino oscillation searches, for example at E776 and E531, as well as the present 



CHORUS ||3T[ and NOMAD as well as the future experiments being discussed at 
present. 



II. ATMOSPHERIC NEUTRINO FLUXES 



A contemporary calculation of the atmospheric neutrino fluxes consists of a Monte 
Carlo procedure that folds the measured energy spectra and chemical composition of the 
cosmic ray flux at the top of the atmosphere with the properties of hadronic interaction 
with the light atmospheric nuclei. Since the properties of the secondary mesons are 
extremely well known, the accuracy of the calculation is determined by the uncertainty of 
the two sets of assumptions - about the primary cosmic ray flux and about the hadronic 
interactions on light nuclei. 

In order to avoid the uncertainty in the absolute magnitude of the cosmic ray flux 
experiments measure the ratio of electron to muon neutrinos, which is very stable in 
different calculations. The absolute normalization of the atmospheric neutrino flux is still 
very important for the interpretation of the observed muon neutrino deficit. If it turns out 
that the measured numbers of electron-like interactions agrees with the predictions and 
there is an absolute deficit of muon-like interactions, the causes must be disappearance. 
If experiments measure the right amount of and an excessively large number of Ve there 
must be a reason for z/g appearance, such as z/^ z/g oscillations or a background process 
that generates Ve or events in the detectors. 



We use the new neutrino flux calculations of the Bartol group |10] which are per 



formed with an updated version of the TARGET interaction event generator. The cosmic 



ray flux model is discussed in detail in Ref. [10|. The treatment of the hadronic collisions 



is very similar to that in earlier calculations There are only minor improvements 

in the treatment of the resonant region for low energy collisions and in the cross section 
for production of KK pairs above 1000 GeV. The improvements in the low energy (2-3 
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GeV) range affects slightly the ffuxes of 100 - 300 MeV neutrinos, while the kaon spectra 
at high energy change the neutrino to anti-neutrino ratios above 100 GeV and thus mostly 
affect the predictions for the flux of upward going neutrino induced muons. 

In the absence of geomagnetic effects the fluxes of GeV neutrinos are practically the 



same as those of Ref. |]13[. Much more significant difference is introduced by the improved 
treatment of the geomagnetic effects [T^. The probability for low rigidity cosmic rays to 
penetrate to the atmosphere and to produce neutrinos is calculated using a realistic model 
of the geomagnetic field and accounting for the shadow of the Earth. As a result the 
neutrino fluxes at experimental locations with high geomagnetic cutoff, such as Kamioka, 
are significantly lower than in Ref. []T^. At high geomagnetic latitude the new fluxes are 
comparable to the original ones. 

This new set of fluxes, as well as the original one, belongs to a group (together with 



the calculation of Ref. |15|) of atmospheric neutrino flux predictions of relatively high 
magnitude. The expected magnitude of the atmospheric neutrinos was discussed by the 
authors of different predictions ||16| who identified the reason for the differences in the 
treatment of the nuclear target effect in the hadronic collisions in the atmosphere. Calcu- 
lations that assume that pion multiplicities in pp and pAir collisions are similar [|1^ predict 
low neutrino flux magnitude. The event generator TARGET produces pion multiplicity 
that is higher by a factor of ~1.6 in pAir interactions above the resonant region. 

The muon fluxes at different atmospheric depths generated with the same code as 
the new neutrino fluxes were compared to the measurements of the MASS experiment ||18|| . 
The predicted altitude profile of muons of energy above 1 GeV agrees with the measured 
one extremely well. 



III. ATMOSPHERIC NEUTRINO CROSS SECTIONS AND EVENT 

DISTRIBUTIONS 

For each experiment the expected number of //-like and e-like events, iV°, a = fi,e, 
in the absence of oscillations can be computed as 

M a , ^a{h,cos9,,E,)-—e{E^)dE,dE^d{cos9,)dh . (1) 
dEyd[cosUy) dEa 

Here Ei, is the neutrino energy and $q, is the flux of atmospheric neutrinos of type a = 
fi, e; Ea is the final charged lepton energy and e{Ea) is the detection efficiency for such 
charged lepton; a is the neutrino-nucleon interaction cross section, u N —>■ N' i; 6,y is the 
angle between the vertical direction and neutrinos (cos 6^=1 corresponds to the downward 
direction). For some experiments, such as Frejus, we also include neutral current events 
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which are misidentified as charged current ones. In eq. ([l|) is the distribution of h 
which is the slant distance from the production point to the sea level for a type neutrinos 
with energy and zenith angle 6^. We took the distribution from ref. [|1^] which is 
normalized as 

J Ka{h, cos 9^, E^)dh = 1. (2) 

As discussed in Sec. II, the neutrino fluxes, in particular in the sub-GeV range, 
depend on the solar activity. In order to take this fact into account we use in eq. (0) and 
also in sec. IV the averaged neutrino flux defined as follows, 

where (J)^**^ and are the atmospheric neutrino fluxes when the sun is most active 
(solar maximum) and quiet (solar minimum), respectively. The coefficients Cmax and Cmin 
(= I — Cmax) are determined according to the running period of each experiment assuming 
that the flux changes linearly with time between solar maximum and minimum. This is 
a first order correction for the solar modulation of the primary cosmic ray fiux which has 
not been included in previous analyses. 



A. Cross Sections 



In order to determine the expected event rates for the various experiments we use 
the neutrino-nucleon cross sections presented in Fig. |l|. We consider separately the 
contributions to the cross section from the exclusive channels of lower multiplicity, quasi- 
elastic scattering and single pion production, and include all additional channels as part 
of the deep inelastic (DIS) cross section |l20| : 

(^^^ = (yQE + cTi^ + (yois ■ (4) 



The quasi-elastic cross section for a neutrino with energy is given by ||21|| : 



da, 



QE 



d\q^ 



vn 



i-p) 



M^Glcos'^e^ 
SirE?, 



u 



M2 



u 



M4 



(5) 



where s — u = AME^ + q^ — mj, M is the proton mass, is the charged lepton mass and 
g2 is the momentum transfer. For up £^n, the same formula applies with the change 
A2 —A2. The functions A\^ A2, and A3 can be written in terms of axial and vector 
form factors: 
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where we have neglected second order currents and we have assumed CVC With this 



assumption all form factors are real and can be written as 

-2 



1 - 



Fa = F^(0) 1 




1 + /Xp - fin) 



AMI 



(7) 



fip and fin are the proton and neutron anomalous magnetic moments and the vector mass, 
My = 0.71 GeV^, is measured with high precision in electron scattering experiments. The 
largest uncertainties in this calculation are associated with the axial form factor. In our 
simulation we use -^4(0) = —1.23 which is known from neutron beta decay. The axial 
mass used by the different collaborations varies in the range M\ = 0.71 — 1.06 GeV^. 

So far we have neglected nuclear effects. The most important of such effects is due 
to the Pauli principle. Following Ref. PT| we include it by using a simple Fermi gas model. 
In this approximation the cross section of a bound nucleon is equal to the cross section 
of a free nucleon multiplied by a factor (1 — N"^D). For neutrons 



D = Z 



for 2z < u — V 



D = Ia\1- f (m2 + v^) + ^+ - v^] for u-v <2z <u + v 



D = 



for 2z > u + V 



with z 



(g2 + m2)2/(4M2 



/{2k 



u 



{2N/AY/^, and v = {2Z/Ay/^. Here 



A,Z,N are the nucleon, proton and neutron numbers and kf is the Fermi momentum, 
kf = 0.225(0.26), for oxygen (iron). For protons, the same formula applies with the 
exchange N Z. The effect of this factor is to decrease the cross section. The decrease 
is larger for smaller neutrino energy. For energies above 1 GeV the nuclear effects lead to 
an 8 % decrease on the quasi-elastic cross section. 



For single pion production we use the model of Fogli and NarduUi p2[ which includes 
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hadronic masses below W = lA GeV. Deep inelastic cross sections are usually described 
in terms of the variables y = 1 — E^/ Ey and x = —q^/{2MEyy). In the parton model 



dxdy 



( \ 



V 



V ) 



Gpsx 

4:71 



F,TF3 + {Ei±Fs){l~yf 



(9) 



where Fi and F^ are given in terms of the parton distributions. For isoscalar targets 
Fi = 2^{qi + (ji) and F3 = ^(gi — qi). In order to avoid double counting we follow 

i i 

the approach of Ref. ||20| and we integrate the deep inelastic contribution in the region 
W >Wc which implies 2MEyy{l - x) > - where Wc = 1.4 GeV. 

The final necessary ingredient are the detector efficiencies given by the experiments. 
These are, in general, functions of the incident neutrino energy and the detected lepton 
energy and flavour. We took these efficiencies from refs. for Kamiokande sub-GeV 
and 1MB, for Frejus, |^ for Nusex. The efficiencies for the Kamiokande multi-GeV 
and e-like events for Soudan2 are provided by the experimentalists and for Soudan2 /i-like 
events the efficiencies are determined in such a way that the energy distributions are well 
reproduced. For the Super-Kamiokande we are making some approximation based on the 
information provided also by the experimentalists, as discussed below. 



B. Event Distributions 



In order to verify the quality of our simulation we compare our predictions for 
the energy distribution of the events with the Monte Carlo simulations of the different 
experiments in absence of oscillation. In Fig. ^we show our predictions superimposed with 
those from the experimental Monte Carlo for the sub-GeV experiments, Kamiokande sub- 
GeV 0, 1MB [|], Frejus 0, Nusex ||] and Soudan II ||. We can see that the agreement 
is very good. No additional normalization of the event rates has been performed. 

Similarly, in Fig. |^ we show the distribution of the fully contained electron-like 
events and fully and partially contained muon events for the Kamiokande multi-GeV 
sample compared with the experimental Monte Carlo prediction given in Fig.2 of Ref 
Some comments are due. In order to obtain these distributions we have used detailed 
experimental efficiencies of Kamiokande for detecting fully contained and partially con- 



tained electron and muon events One must take into account that the Monte Carlo 
distributions given in Fig.2 of Ref were generated using the fluxes of Honda et al. [|15 



while we are using the fluxes of Gaisser et al. []T3[ . Thus we have an absolute prediction for 
the number of events for Kamiokande multi-GeV data and for their energy distribution 
which is obtained under the same assumptions for the cross sections and neutrino fluxes 
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than the other sub-GeV experiments. For the sake of comparison we also show in Fig. ^ 
(upper two panels) the angular distribution of the events for Kamiokande multi-GeV in 
the absence of oscillations as obtained from our calculation. 

In Fig. ^ we plot as in Fig. ^ the expected energy distribution also for Super- 
Kamiokande sub-GeV and multi-GeV data. We also plot in Fig. ^ the angular distribu- 
tion for Super-Kamiokande sub-GeV (middle two panels) and multi-GeV data (lower two 
panels). We have used, as an approximation, the preliminary acceptances of Super- 
Kamiokande for 325.8 days for fully contained, 293 days for partially contained events as 
detection efficiencies for final leptons, for sub-GeV as well as multi-GeV data. In order 
to obtain the angular distribution of expected events for in the multi-GeV range we have 
assumed that the lepton direction is the same as the incident neutrino direction. Actually 
for the Kamiokande multi-GeV data, the average angle between the incident neutrino and 
the lepton direction is about 15°. In our calculation we have simulated this difference by 
smearing the angular distribution with a Gaussian distribution with one-sigma width of 
15°. As seen in Fig. ^ the effect of this approximation is small. At this point it is worth 
noting that the angular distribution for multi-GeV electrons in the Super-Kamiokande 
sample is flatter than in the Kamiokande data. The main reason for this zenith-angle 
shape difference is due to the smaller selection efficiency for the 1-ring e-like events at 
high energy in the super-Kamiokande analysis [|23| . As a result the mean neutrino energy 



was shifted to lower value and the mean angle between the incident neutrino and the lep- 
ton direction becomes larger. We have simulated this effect by increasing the one-sigma 
width of the smearing Gaussian to 25° for the super-Kamiokande multi-GeV electrons 
which effectively flattens the angular distribution as seen in Fig. |[ 

On the other hand, for events in the sub-GeV range we have carefully taken into 
account the difference between the incoming neutrino angle and the detected charged 
lepton scattering angle which is a function of the incoming neutrino energy. As can be 
seen in Fig. ^ this leads to a much flatter expected angular distribution for the sub-GeV 
neutrinos, in agreement with the prediction from the experimental MC. 



We also estimate the expected the ratio in the absence of oscillation as 

-)0 ^ fV 



where A*"^ and are computed by eq. (|T]). In Table || we present our prediction for the 
expected ratio in the absence of oscillations for the various experiments and compare it 
with the expected MC results [1]-^]. Table |I| also displays our prediction for the expected 
ratio for the Kamiokande multi-GeV and Super-Kamiokande zenith angle distribution. 
We see that the agreement between Rff/^ and our prediction are very good for most 
of the experiments. 
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IV. ATMOSPHERIC NEUTRINO DATA FITS 



We now consider the simplest interpretation of the atmospheric neutrino anomaly 
in terms of the neutrino oscillation hypothesis. For definiteness we assume a two-flavor 
oscillation scenario, in which the oscillates into another flavour either z/^ —>■ or 

I/^ Vr- 



A. Data Analysis Procedure 



In the presence of two-flavour neutrino oscillations, the expected number of (i and 
e-like events, Na, a = /x, e is given by 

N, = N^iP,,) + Nl (P,^) , iVe = Ni{P,,) + iVj,(P^,) , (11) 

where 

^°/3= / ,^ " ^ . Kp{h,cos9,,E,)-—e{Ep)dE,dEpd{cos9,)dh (12) 

and 

{Pap) = t4- / jr^t^" n . i^i3{h,cose^,,E^)Pap-^e{Eii)dE^dEiid{cose^)dh . (13) 

Here Pap is the oscillation probability of — > z/^ for given values of E^^ , cos 6^ and /i, 
i.e., P^p = P{vp z/c,; S^.^, cos6'j„ h) 

We note that for the z/^ — > channel. Earth matter effects lead to oscillation 
probabilities which are different for neutrinos and anti-neutrinos. Therefore, we separately 
compute Pafi in eq. ([T3| ) in order to correctly estimate the expected number of events in 
each experiment (see sec. IV. B for more discussion). 

When combining the results of the experiments we do not make use of the double 
ratio, i?^/e/-R^f , but instead we treat the e and /x-like data separately, taking into account 
carefully the correlation of errors. Following ref. we define as 



= E(^^f - Nf-') ■ (^L, + a,\_^)7] ■ {Nr" - Nf-y), (14) 



theory J 
I,J 



where / and J stand for any combination of experimental data set and type of events 
considered, i.e, I = {A, a) and J = (P,/3) where, A,B = Frejus, Kamiokande sub-GeV, 
1MB,... and a,P = e,fi. In eq. (|14D ]\f^^^'^°^y is the predicted number of events calculated 
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by eq. (0) whereas Nf°'*°' is the number of observed events. In eq. (|l^) a^^ta crli^eory 



are the error matrices containing the experimental errors and the MC errors respectively. 
They can be written as 

a]j = aM)Pc.p{A.B)cyp{B), (15) 

where pajsiA, B) stands for the correlation between the a-like events in the A-type experi- 
ment and /3-like events in S-type experiment, whereas (Ta{A) and crp{B) are the errors for 
the number of a and /5-like events in A and B experiments, respectively. The dimension 
of the error matrix varies depending on the combination of experiments included in the 
analysis. For each individual experiment, the error matrix has dimension 2x2 whereas 
for the full experimental data set with binning (20 data for each flavor) its dimension is 
40 X 40. 

With this procedure of treating separately the e-like and /i-like data with correlation 
of errors, we avoid the non-Gaussian nature of double ratio, as pointed out in ref. 



We compute pap{A, B) as in ref. |2^. A detailed discussion of the errors and cor- 
relations used in our analysis can be found in the Appendix. In Table |T| we show the 
values of and the confidence level in the absence of oscillation. In our analysis, we have 
conservatively assumed 30% uncertainty regarding to the absolute neutrino flux, in order 
to generously account for the spread of neutrino flux predictions in different calculations 

n 



Next we minimize the x function in eq. ([T^ ) and determine the allowed region in 



the sin^ 29 — /S.rn? plane, for a given confidence level, defined as, 

x' = xLn + 4.61 (9.21) for 90 (99)% C.L. (16) 



B. ^ Uf. channel 

The results of our fit of atmospheric neutrino data obtained at the various indi- 
vidual water- Cerenkov and iron calorimeter detectors for the — * z/g channel are shown 
in Fig. |[ The left-pointing arrows in the left figure in Fig. ^ correspond to the nega- 
tive results of Frejus and Nusex. So far we have not included in the above analysis the 
constraints that arise from the inclusion of the angular dependence of the data in the 
Kamiokande multi-GeV data as well as the Super-Kamiokande data. In the right hand 



^For a brief discussion of the effect of the assumed flux uncertainties, see section V 
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panel of Fig. ^ we show how the binned results of Kamiokande and Super-Kamiokande 
give rise to a region of oscillation parameters that cuts out the large Am^ values. More- 
over one can see that the Super-Kamiokande binned sub-GeV data yield a somewhat lower 
value of Am^ than the multi-GeV data. 

The effects of combining all atmospheric neutrino data from the various experiments 
for the z/^j channel are shown in Fig. |^. This figure show the allowed v^. 

oscillation parameters for all experiments combined at 90 and 99 % CL. For comparison 
we have also plotted in Fig. ^ the presently excluded region from reactor experiments, 



Krasnoyarsk |28|, Bugey [59] and the recent Chooz long-baseline result [12 



We have so far neglected Earth matter effects [^, both in Fig. |]and Fig. 0. In order 



to take into account the earth matter effect in our analysis we have separately computed, 
by numerical integration, the oscillation probabilities, P[y^i —>■ Ve) = P{^e ^ ^'^t) and 
P{i^fi ^e) = P{i^e — ^ i^fi)- This is necessary, since the matter effect distinguishes 
neutrinos from anti-neutrinos. We have used the approximate analytic expression for the 
electron density profile in the Earth obtained in ref. [0. In order to save computation 



(CPU) time we have neglected the matter effect for neutrino oscillation parameters in the 
range, 

> 10-^^eV, (17) 



E 

since the maximum value of the matter potential (at the Earth center) is at most, 

Vmatter ~ IQ-^^eV, (1^ 



and the matter effect on the probability is small if condition eq. ([T7|) is satisfied. 

In Fig. H we show the allowed z/^ z/g oscillation parameters for each individual 
experiment including Earth matter effects. Unlike the previous case where matter effects 
were neglected, a noticeable new feature in this case is that the Super-Kamiokande multi- 
GeV data now allows large Am^ values, even if binning is taken into account. The allowed 
i^ti oscillation parameters for the Super-Kamiokande binned data combined at 90 and 
99 % CL including Earth matter effects in shown in Fig. An interesting feature to note 
here is that by adding the matter effects the allowed regions lie higher in Am? than when 
matter effects are neglected. This is because for smaller Am^, i.e. when Am? cos26/2E 
is much smaller than Vmatter the effective conversion amplitude sin^ 29m where 9m is the 
mixing angle in matter, is smaller than that of the vacuum one, i.e. sin^ 29. In other 
words, in this region matter suppresses the conversion and it becomes harder to fit for 
the atmospheric neutrino anomaly. 

The allowed z/^ Vf, oscillation parameters for all experiments combined at 90 and 
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99 % CL including Earth matter effects is shown in Fig. |T0|. Again one can see that by 
adding the matter effects the allowed regions lift higher in Am^ than when matter effects 
are neglected. We found the best fit point at (sin^ 26, Am^) ~ (0.97, 2.6 x 10~^ eV^) where 
Xmin — 62.7 for 40 degrees of freedom. We would like to point out that at this stage the 
weight of the experiments with negative results (NUSEX and Frejus) is small enough not 
to modify the per degree of freedom (Xmin — ^^.2 for 36 degrees of freedom when these 
experiments are removed). 

It is instructive at this stage to compare the region determined by the atmospheric 
neutrino data fit with the presently excluded region from reactor experiments p9|. The 
inclusion of the matter effects becomes especially relevant when one compares with the 



long baseline reactor neutrino data, such as the recent data of Chooz |T^. One sees that 
at 90% CL the i/^ to z/g oscillation channel is ruled out as a solution of the atmospheric 
neutrino anomaly. 



C. ^ Vt channel 



The results of our x fit of atmospheric neutrino data obtained from the data of 



individual experiments for the channel are shown in Fig. [IT]. The allowed 

regions for each experiment lie to the right of the corresponding labelled line, except for 
the negative Frejus and Nusex experiments. In the left figure in Fig. we have not 
included the constraints that arise from the inclusion of the angular dependence of the 
data in the Kamiokande multi-GeV data as well as in the Super-Kamiokande data. 

It is instructive to compare the results obtained for the Kamiokande data with 
those obtained by including the recent Super-Kamiokande data. In Fig. |12| we show the 
allowed z/^ — ^ oscillation parameters for Kamiokande and Kamiokande plus Super- 
Kamiokande combined. Some features are worth remarking. For example, the inclusion 
of the unbinned Super-Kamiokande data to the corresponding Kamiokande data leads to 
the exclusion of large mixing in the large Am^ region. On the other hand the inclusion 
of Super-Kamiokande binned data leads to a substantially smaller region that obtained 
from the Kamiokande full data sample, reflecting a real improvement. 



In Fig. O we give the allowed z/^ — Vr oscillation parameters for Super-Kamiokande 
combined at 90 and 99 % CL, while in Fig. [1^ we display the allowed z/^ oscillation 



parameters for all experiments combined at 90 and 99 % CL. By comparing Fig. O with 
Fig. |14| one can see the weight of the Super-Kamiokande data sample in the total data 
sample collected by all experiments. We find that the best fit points lie at (sin^ 2^, Am^) ~ 
(1, 1.3 X lO^^eV^) with Xmm = 14.4 for the 20 degrees of freedom, for Super-Kamiokande 
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only and {sin^ 26 , Arn^) ~ (1,1-2 x lO^^eV^) with xLin = 66.6 for the 40 degrees of 
freedom, for all combined. The global fit to all experiments is still slightly better for the 
channel. However, the difference between the quality of the fit for both channels 
is smaller now than in the pre-Super-Kamiokande era, due to the angular distribution of 
Super-Kamiokande multi-GeV which strongly favours the channel. 

The result of including the information on the zenith angle distribution of the events 
in the v^^ to v.^ fit is clearly to cut the large values of Am^, as can be seen in all figures, 
namely Fig. ^ Fig. |g. Fig. |T3| and Fig. ^. 

One point worth noting is that the inclusion of the new Super-Kamiokande data 
produces a downwards shift in the (sin^ 2^, Am^) region, when compared with pre-Super- 
Kamiokande fits. The importance of the information obtained from the analysis of the 
atmospheric neutrino data analysis in relation to the results from accelerator experiments 



such as E776 and E531 and CDHSW |g as well as CHORUS plus NOMAD combined 
limits and the prospects for the future experiments being discussed at present can 
be appreciated in Fig. 0. One sees that the long-baseline experiments planned at KEK 
(K2K) H, Fermilab (MINOS) H and CERN (NOE |6l and ICARUS ^) fall short 
in sensitivity to probe the z/^ to Vr oscillation parameters. This is in contrast with the 
situation in the pre-Super-Kamiokande days. From this point of view experiments such 
as ICARUS and a re-design of experiments such as MINOS would be desired in order to 
enhance their sensitivity in testing the atmospheric neutrino anomaly. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have considered the impact of recent experimental results on at- 
mospheric neutrinos from Super-Kamiokande and Soudan2 as well as recent theoretical 
improvements in fiux calculations and neutrino-nucleon cross sections on the determina- 
tions of atmospheric neutrino oscillation parameters, both for the — > and z/^ — > 
channels. The new Super-Kamiokande data cause a downwards shift in the (sin^ 2^, Am^) 
region, when compared with pre-Super-Kamiokande results. We have also compared the 
results obtained in our fits of atmospheric neutrino data with previous results, as well 
as with the constraints following from laboratory searches for neutrino oscillations, both 
at accelerators and reactors. For example we have seen that the v^. oscillation 

hypothesis is barely consistent with the recent negative result of the Chooz reactor |T^]. 
The sensitivity attained in atmospheric neutrino observations in the channel is 

also compared with those of accelerator neutrino oscillation searches, for example at the 
present CHORUS and NOMAD as well as at the future experiments being discussed at 
present. Specially interesting from our point of view are the long-baseline experiments 
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planned at KEK (K2K), Fermilab (MINOS) and CERN (NOE, ICARUS). However, due 
to the lowering of the allowed (sin^ 29, Am^) region, it is not clear whether a re-design 
is needed in some of these experiments, for example MINOS, in order to enhance their 
sensitivity in testing the atmospheric neutrino anomaly. 

Note that, throughout this work we have assumed a rather generous error in the 
absolute fluxes of atmospheric neutrinos. We have investigated to some extent the effect 
that a reduced error in the fluxes would have in the determination of neutrino oscillation 
parameters from the present atmospheric neutrino data. Wc have found no significant 
effect in the shape of allowed region when we changed the assumed error in the fiuxes from 
30 % to 20 %. However, wc have noticed a somewhat significant effect of a more accurate 
ratio of muon-to-electron events. We have found, for example, for Super-Kamiokande, 
that when we decrease the error in the muon-to-electron-type event ratio from 5 % (10 
%) for unbinned (binned) data to 3 % (6%) the allowed region shrinks by about 10 to 15 
% sin^ 29 close to 0.7 or so. There is hardly any effect in the Am^ range determination. 
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APPENDIX: CORRELATION OF ERRORS 



Here we present the errors and correlations used in our analysis. In Table |^ we 
display the errors and correlations p^e(^, ^) for all the experiments. Data errors and 
correlations contain the experimental statistical errors as well as those due to misidenti- 
fication as quoted by the experiments. In order to compute (j^^^°^y we take into account 



25| , the flux uncertainty, the MC statistical errors (which depend on the number of sim- 
ulated MC events) as well as the cross sections uncertainty. The flux uncertainty is taken 
to be 30% whereas MC statistical errors are estimated under the assumption that the /x 
and e-like events follow a binomial distribution. Nuclear cross section uncertainties are 
taken to be 10 % for all the experiments except for Soudan2 for when we used the values, 
7.5% and 6.4% for e-like and /x-like events, respectively 0. 

Data errors between different experiments are assumed to be uncorrelated, 
pit {A, A) = 1 (« = e,/i) for all A, 

while the theory correlations between different experiments (i.e., for A ^ B) are obtained 
as follows, 

flux flux 

theory, . ry^ _ flux a f3 \f ry =^ H 

Paf3 1^5 -DJ - Pa/3 ^ ^theory ^j^^ ^theory ^j^^ U a ^ p, 

flux flux 

theory f A td\ _ flux a P ;f ^ _ O 

Pq/3 1^5-DJ-Pq/3 ^ theory,.-. theory f rj, UO-p, 

where = a/;'"^' = 30 % and pi^^"" = 1.000 for a = p. For a =^ p, we use p{jp^ = 0.986 
(0.944) as determined from the relation 

H'P' = K'n' + - 2p,ei^i'ni^['n 

after imposing that the uncertainty in the flavor ratio, crj^'l^f = 5%(10%) for unbinned 
(binned) case Furthermore we assume that there is no correlation between the sub- 



GeV and multi-GeV data. 

We note that for both sub-GeV and multi-GeV data, in general, p*^p"^^{A, B) is not 
symmetric under the exchange of the flavor labels a and (3 or the experimental labels , 
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B) ^ p'^r^iA, B) if a ^ /3, 
but it is symmetric under simultaneous exchange of both kinds of labels a, (3 and A, B: 

Pari^^B)^p'^ry{B,A). 
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TABLES 



Experiment 




Super-Kamiokande sub-GeV 


0.635 ± 0.035 ± 0.053 


Super-Kamiokande multi-GeV 


0.604 ± 0.065 ± 0.065 


Soudan2 


0.61 ±0.14 ±0.07 


1MB 


0.55 ±0.11 


Kamiokande sub-GeV 


0.6 ± 0.09 


Kamiokande multi-GeV 


0.59 ±0.1 


Prejus 


1.06 ±0.23 


Nusex 


0.96 ± 0.3 



TABLE I. Results from the atmospheric neutrino experiments 
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^MC 


V 

IvF 


Prejus 


1.9 


1.8 


Super-Kamiokande (sub-GeV) 


1.6 


1.6 


Kamiokande (sub-GeV) 


1.55 


1.6 


Binl 


1.7 


1.6 


1MB 


1.1 


1.1 


Bin2 


1.6 


1.5 


Soudan2 


1.05 


1.1 


Bin3 


1.5 


1.5 


Nusex 


1.9 


1.8 


Bin4 


1.5 


1.6 


Kamiokande (multi-GeV) 


2.3 


2.4 


Bin5 


1.7 


1.5 


Binl 


3.1 


3.1 


Super-Kamiokande (multi-GeV) 


3.2 


3.0 


Bin2 


2.4 


2.4 


Binl 


3.8 


3.4 


Bin3 


2.1 


2.0 


Bin2 


2.8 


2.8 


Binl 


2.4 


2.4 


Bin3 


3.2 


2.8 


Bin5 


3.2 


3.2 


Bin4 


2.9 


2.8 








Bin5 


4.2 


3.5 



TABLE II. Our predictions for the ratio {N^/N^) in the absence of oscillations compared 
to the MC expectations {Njf-^ /N^'-') from each experimental group. 
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Experiment 




C. L. (%) 


Prejus 


0.56 


24.4 


1MB 


8.4 


98.5 


Soudan2 


5.7 


94.2 


Nusex 


0.39 


17.7 


Kamiokande sub-GeV 


12.5 


99.8 


Kamiokande multi-GeV unbinned 


8.7 


98.7 


Kamiokande multi-GeV binned 


18.2 


94.8 


Super-Kamiokande sub-GeV unbinned 


21.5 


99.7 


Super-Kamiokande sub-GeV binned 


27.2 


100.0 


Super-Kamiokande multi-GeV unbinned 


10. 


99.3 


Super-Kamiokande multi-GeV binned 


27.9 


99.8 



TABLE III. Values of and confidence level for each experiment in the absence of oscilla- 
tions. For unbinned data the number of degrees of freedom is 2 while for combined binned data 
is 10. 
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HjXpcIlIlicIll I J 




(7g (/oj /? 


data (A A\ 
ixe K-^i 


^theory 


^theory ptf/o^2/(A, A) 


irejus 


10.5 


17.9 


—0.021 


31.7 


31.9 


0.951 


Kam sub-GeV 


7.1 


7.0 


-0.081 


31.7 


31.8 


0.975 


1MB 


8.9 


7.5 


-0.374 


36.1 


36.1 


0.947 


Nusex 


18.4 


27.2 


-0.050 


31.7 


31.9 


0.950 


boudani 


13. 5 


11.0 


— O.loo 


OA O 

30.8 


31.1 


0.960 


buper-Kam sub-GeV (unbmnedj 


4.9 


4.9 


—0.042 


31.6 


31.7 


0.978 


buper-Kam sub-GeV bml 


9.4 


8.3 


—0.013 


31.7 


31.8 


0.936 


buper-Kam sub-GeV bm2 


9.L) 


9.4 


—0.012 


31.7 


31.8 


0.935 


buper-Kam sub-GeV bm3 


9.0 


8.4 


—0.013 


31.7 


31.8 


0.936 


Super-Kam sub-GeV bin4 


O.D 


O 1 

y.i 


— U.Uio 


t5i. / 


31.8 


0.936 


^ 1 1 OT '1 iT"i c;i 1 l"i f-loA ^ In 11 
O 1 1 1 J I I -t\ f 1 1 ii n 1 1 U v Tl V V )\ \A. - ) 


o . v> 






■^1 7 


31.8 


0.935 


Kam muiti-GeV (unbmnedj 


9.D 


11.0 


A AO O 

—0.038 


31.7 


32.0 


0.965 


Kam muiti-GeV bml 


24.0 


22.2 


—0.008 


31.8 


33.9 


0.840 


Kam multi-GeV bm2 


22.8 


23.3 


—0.008 


31.9 


32.9 


0.869 


Kam multi-GeV bin3 


18.2 


19.1 


-0.012 


31.8 


32.6 


0.889 


Kam muiti-GeV bm4 


18.8 


22.8 


-0.009 


31.9 


32.8 


0.878 


Kam multi-GeV bin5 


17. z 


66.0 


— U.iJUT 


oi.y 


33.8 


0.838 


buper-Kam muiti-GeV (unbmnedj 


7.0 


9.5 


— U.05d 


31.0 


31.9 


0.972 


buper-Kam muiti-GeV bml 


17.9 


21.9 


—0.010 


31.7 


32.4 


0.911 


buper-Kam multi-GeV bm2 


17.4 


18.1 


—0.013 


31.7 


32.1 


0.920 


Supcr-Kam multi-GeV bin3 


12.4 


18.9 


-0.017 


31.7 


32.1 


0.923 


Super-Kam multi-GeV bin4 


12.2 


17.1 


-0.019 


31.7 


32.1 


0.919 


Super-Kam multi-GeV bin5 


12.7 


19.8 


-0.016 


31.7 


32.4 


0.909 



TABLE IV. Errors and correlation for both observed data and theory (MC) samples. 
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FIG. 1. Neutrino-nucleon cross section used in this paper. 
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FIG. 2. Expected energy distribution of Sub-GeV events (histogram) compared with our 
prediction (full line). 
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FIG. 3. Expected neutrino energy distribution of Kamiokande Multi-GeV events (dashed 
histogram) compared with our prediction (fuh histogram) . 
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prediction. Both our prediction and the MC prediction are based on the same flux calculations 

ic 



26 



5G 



(D 40 

> 

(D 

-7 30 



20 
10 




II 

- Kam 


mult' 


-deV 




1 


-FC e 




+ 


















-I- 
















T - 


" 1 1 1 1 


1 


1 1 


1 1 


III" 


1 


G.5 





0.5 





m 
c 

(D 
> 
0) 



cos(0^ 



20 



Kam Mulii-GeV 
FC + PC /J. 



-1 



-0.5 



0.5 1 

cos(O) 



5^200 



100 



SKam Sub — GeV 
1 —ring e 



- tn 



CD 

> 200 



100 



SKam Sub — GeV 
1 —ring /i, 



to 



60 



5j 40 



20 



-0.5 

^ I I I I I I I I I I 
- SKam Multi-GeV 

_ 1 —ring e 
I 



0.5 

cos(0^ 



> 
CD 



mi 



50 



100 



50 



-0.5 



-SKam Multi-GeV 
1 — ring + PG fi 



0.5 

cosfO 







— I — , 




^ 





-1 -0.5 



0.5 1 

cos(O) 



-1 -0.5 



0.5 1 

cos(O) 



FIG. 5. Expected angular distribution of Kamiokande multi-GeV events and Su- 
per-Kamiokande events (dashed histogram) obtained by Monte Carlo simulation by the ex- 
perimental group compared with our predictions (full histogram) and the experimental data. 
We note that in these figures the MC prediction is based on Honda et al. fluxes whereas 



ours is based on Bartol fluxes |fO| normalized to the total number of expected events with the 
Honda MC. 
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FIG. 6. Allowed — > oscillation parameters at 90 % CL for each individual experiment 
neglecting Earth matter effects. 
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FIG. 7. Allowed oscillation parameters for all experiments combined at 90 and 99 % 

CL neglecting Earth matter effects. For comparison we also plot the presently excluded region 
from reactor experiments. 
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FIG. 8. Allowed Ue oscillation parameters at 90% CL for each individual experiment 

including Earth matter effects. 
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FIG. 9. Allowed i/g oscillation parameters for Superkamiokande experiment combined 

at 90 and 99 % GL including Earth matter effects. 
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FIG. 10. Allowed v^, oscillation parameters for all experiments combined at 90 (solid) 

and 99 % CL (dashed) including Earth matter effects. For comparison we also plot the presently 
excluded region from reactor experiments. The cross represents the best fit point. 
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FIG. 11. Allowed oscillation parameters at 90% CL for each individual experiment. 
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FIG. 12. Allowed Ut oscillation parameters at 90% CL for Kamiokande and 

Kamiokande plus Super kamiokande combined. 
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FIG. 13. Allowed — > Vt oscillation parameters for Superkamiokande combined at 90 and 
99 % CL. 
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FIG. 14. Allowed v^j, Vr oscillation parameters for all experiments combined at 90 and 99 

% CL. For comparison we also display the presently excluded region from accelerator experiments 
CDHSW and CHORUS+NOMAD and future long baseline experiments. 
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